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In this paper, we study some interesting properties of a spherically symmetric oscillating soliton 
star made of a real time-dependent scalar field which is called an oscillaton. The known final 
configuration of an oscillaton consists of a stationary stage in which the scalar field and the metric 
coefficients oscillate in time if the scalar potential is quadratic. The differential equations that arise 
in the simplest approximation, that of coherent scalar oscillations, are presented for a quadratic 
scalar potential. This allows us to take a closer look at the interesting properties of these oscillating 
objects. The leading terms of the solutions considering a quartic and a cosh scalar potentials 
are worked in the so called stationary limit procedure. This procedure reveals the form in which 
oscillatons and boson stars may be related and useful information about oscillatons is obtained 
from the known results of boson stars. Oscillatons could compete with boson stars as interesting 
astrophysical objects, since they would be predicted by scalar field dark matter models. 

PACS numbers; 95.30.Sf, 98.80.-k 



I. INTRODUCTION 

Being the nature of dark matter one of the most in- 
triguing open questions in physics nowadays jl]-^, it is 
not strange that the number of proposals trying to solve 
it by considering exotic matter has increased in recent 
years. Among them, models of bosonic dark matter seem 
promising and have been widely studied, both at cosmo- 
logical and galactic scales p|-p8[|. On one hand, complex 
scalar field boson stars are the most studied as objects of 
astrophysical interest, see for instance poj-^ and many 
references there in. On the other hand, bounded objects 
formed by real scalar fields has not been studied that 
deeply yet [|||jl|jl|,|l|,p,|4|,|9). 

The study of gravitationally bound objects with real 
scalar fields is a very important issue for the scalar field 
dark matter hypothesis at galactic scales ]T^"[l^, but 
most of the recent treatments have been static, both 
Newtonian and relativistic ^ ||,|ll|,|l^ . It is known that 
the only fully relativistic bound object made of real 
scalar field, the oscillaton, is time-dependent in nature 
|p|, P^|JT^j2^ , in opposition to the static treatments that 
can be made for boson stars, i.e., for complex scalar fields. 
Besides the possible interest for cosmologists, from the 
point of view of general relativity, oscillatons are also im- 
portant since they are time dependent, non-topological, 
non singular and asymptotically flat solutions to the cou- 
pled Einstein-Klein-Gordon equations. Their intrinsic 
properties are quite different than those of the known 
static, or even stationary, solutions to the Einstein equa- 
tions. Moreover, their stability and possible formation by 
means of a Jeans mechanism has been studied partially 



before |jT||2|], but more studies are needed in order to 
establish their possible astrophysical role. 

In this paper, our main aim is to make a generic study 
of the solutions to the spherically symmetric coupled 
Einstein-Klein-Gordon equations for a time-dependent 
real scalar field <& endowed with a quadratic potential, 
using a simple approximation. Within this approxima- 
tion, we can manipulate the differential equations, solve 
them numerically and obtain, easily and clearly, some of 
the interesting properties of oscillatons. This elucidates 
the results presented in previous papers about these ob- 
jects. Also, we show how we can extend this study to 
other scalar potentials, which have not been treated be- 
fore but that have been taken into account in scalar field 
dark matter models. 

A brief summary of the paper is as follows. In section 
II, we show how to obtain the equations for an oscilla- 
ton with a quadratic scalar potential, under the assump- 
tion of coherent scalar field oscillations with fundamen- 
tal frequency lu, that is, the scalar field is of the form 
<I>(t, r) = (7(r) cos(wt). While the general form of the 
metric coefficients can be known within this assumption, 
an approximation in Fourier series is done in order to deal 
with the full differential equations. It is discussed in turn 
whether this approximation is enough, even though the 
Klein-Gordon equation for the scalar field has a time- 
dependent damping term. Taking appropriate dimen- 
sionless quantities and boundary conditions, we find some 
non singular and asymptotically flat 0-node solutions. 
Since the solutions are also asymptotically static, in the 
sense that the oscillations of the oscillaton are confined 
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to a finite spatial region, the mass of the configurations 
is the Schwarzschild mass that an observer measures at 
infinity. This ahows us to calculate the critical mass of 
oscillatons. 

In section III, the similarity between a boson star and 
an oscillaton is discussed under the so called stationary 
limit procedure. This limit can be seen as the lowest or- 
der of approximation to an oscillaton in which the differ- 
ential equations become boson star-like. Some numerical 
results are presented for comparison between an oscilla- 
ton and a boson star. The similarity is better for weaker 
fields and shows that oscillatons can become interesting 
astrophysical objects at the same level of boson stars. 

The case of a scalar field with a quartic self-interaction 
is worked within the stationary limit procedure in sec- 
tion IV (A). The resulting equations can also be written 
in a boson star-like form, which permits us to easily ob- 
tain information for oscillatons using the known results 
of boson stars. This limit is also applied to the case of 
a (non-polynomial) cosh-potential in section IV (B), but 
the resulting equations are not boson star-like. However, 
we can use the similarity between cosh and quartic os- 
cillatons to obtain useful information. In section V, we 
summarize the results and give points to be investigated 
in future research. The numerical method used in this pa- 
per is the fourth-order Runge-Kutta method of MAPLE. 



II. OSCILLATONS 

We start by taking the case of a quadratic scalar inter- 
action with spherical symmetry, that was analyzed in . 
The most general spherically-symmetric metric is written 
as 



9af3 



dx°'dx'^ = -e'^^i'dt^ + e^+f'dr^ + r^dn^ (1) 



where v — u[t^r)^ fi — M(ij^) ^^re yet arbitrary functions 
(the units are such c = 1). The energy- momentum tensor 
of a real scalar field ^{t,r) endowed with a scalar field 
potential V{^) is defined as @ 



whose non-vanishing components are 
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and also T^^ — T^2- Overdots denote d/dt and primes 
denote d/dr. These different components are identified 
as the energy density the momentum density 7^$, the 
radial pressure Pr and the angular pressure p± . 



The Einstein equations, Gaf3 = kqTq^, can be written 
as differential equations for the functions v, fi, then 

(i^ + /i)- = Kof^*', (7) 
„' = ^ (e2^$2 + , (8) 

m' = J[l + e"+''('SorV-l)], (9) 

where kq = SttG — Stt/toI,; is the inverse of the reduced 
Planck mass squared. The conservation equations for the 
scalar field energy-momentum tensor (0) are written as 



dV 



= 0, 
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where □ — g'^^Va^ p is the d'Alambertian operator. 
Therefore, we obtain the Klein-Gordon (KG) equation 
for the scalar field <&, 
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Observe that the metric function i' is related to the 
scalar kinetic energy while /i is related to the scalar po- 
tential energy VYp). Actually, eqs. ([7[]ll|) are easier than 
the usual ones g]. Notice that the left-hand side of the 
KG equation has the usual form of the static case and 
then the right-hand side looks like a source term. From 
this, it is not surprising that one can get non-singular 
solutions. 

The quadratic scalar potential is written as V{^) — 
where m is the boson mass. If we choose 
$(i, r) = a{r)(f){t), eq. (O) can be rewritten as 
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This equation is almost separable. The terms on the 
right-hand side suggest that the scalar field oscillates har- 
monically in time with a damping term related to p. Fol- 
lowing the work p|, we just consider that 



(13) 



Ko^{t,r) = 2cr(r)cos (cot), 



where uj is the frequency of the scalar oscillations. 

For simplicity, only this first approximation will be 
used along this paper, in order to analyze and reveal the 
main properties of oscillatons, as the procedure used in 
was somewhat obscure. Besides, this is the most sim- 
ple approximation since we are not including a damping 
term. I will discuss this at the end of this section. 

Eq. M) can be formally integrated up to 



' + /i = (ly + fi)o + raa' cos (2a;i), 



(14) 



with + /i)o an arbitrary function of the r— coordinate 
only. The metric functions fi should then be expanded 
as 
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r) = iyo{f) + cos (2wt), 
/i(t, r) = Mol?") + Mil?") cos (2w<). 



(15) 



Taking into account the Fourier expansion 

oo 

e^^^^™^^''^ =/o(/(a;)) + 2^/„(/(x))cos(2n^?), (16) 

n=l 

where /„(z) are the modified Bessel functions of the first 
kind, we can expand the metric coefficients as 
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In making the expansions (19-2^, the neglected terms 
on the right hand side were those containing cos {^ujt) , 
while the neglected ones in the KG equation were those 
with cos(3w<), cos(5tji). This suggests that the metric 
coefficients should be expanded with even Fourier terms, 
and that the scalar field expansion involves only odd 
Fourier terms. Then, the expansions used in |^ are well 
justified. 

The solution is completely determined by solving eqs. 
(|l9y2^ ). But first, observe that eq. (0) gives the exact 
algebraic relation 



(24) 



The metric coefRcients also oscillate in time but with 
even- multiples of the fundamental frequency lj. 



We will take this last equation to calculate function vi 
.ii^itife^d of integrating eq. (|l^), and then we are left with 
one algebraic and four ordinary differential equations. 



B. Boundary conditions 



A. Differential equations 

For numerical purposes, we perform a change of vari- 
ables as in the boson star case 121^j22l 



X — mr, 
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from which the metric coefRcients are given now by 
grr — e''+^ and gtt — —^^e^~^. Observe that the bo- 
son mass sets both the scale of time and distance. 

The differential equations for the metric functions ( p^ ) 
and the scalar field ( |l3|) appear from equations (|^pT|) by 
setting each Fourier component to zero. In contrast to 
the truncated expansion in B, we can use as much terms 
as necessary in expansions (|17|), i.e., there is no a priori 
truncation. The equations to be solved are 



In order to have nonsingular solutions, regularity at 
a; = requires cr'(O) — and 1^(^,0) + ^J.{t,Q) — 0, so 
j^o(O) = — /io(0) and Mi(0) = — i^i(O). But the latter con- 
dition is satisfied straightforwardly from eq. (p4|). 

Asymptotic flatness require that ct(oo) = /ii(oo) = 0, 
which directly implies i^i(oo) = 0. However, /xo(oo) — 
—vq{oo) ^ because of the change of variables in 

(p8|), and thus exp (z/ — /i)(oo) — giving the 

value of the fundamental frequency w [ElU22], but still 
exp (z/ -|- n){oo) = 1. This also has the advantage that 
for each value of cr(0), we only have two degrees of free- 
dom and then we need only to adjust the central values 
/io(0),/ii(0) to obtain different n-node solutions. How- 
ever, we will deal with 0-node solutions only. 

The central value /ii(0) deserves special attention. 
From eq. (22), we see that its radial derivative is al- 
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ways positive. Thus, the asymptotically flat condition is 
only reached if one chooses Mi(0) < 0. Also, the condition 
I Ml I < 1 is needed for the solutions of eqs. ( p0[p3| ) to con- 
verge, because we can neglect higher terms in expansion 
([TtI). As we will find later (see Fig. the central values 
(t(0), Mo(0), Mi(0) E^re of the same order of magnitude. 

Needless to say, there is only one solution satisfying the 
above boundary conditions for each central value cr(0). 
The procedure to solve numerically (pO|-p4|) is similar as 
for boson stars. 



C. Numerical results 
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Typical metric coefRcients for a 0-node solution with 
a central value (t(0) = 0.2/-\/2 are shown in Fig. |l|. We 
can identify the graphs in the figure by comparing the 
(23) expansion 



where now primes denote d/ dx. 
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with the corresponding expansion in ( p^ for each metric 
coefficient. The general behavior of the metric coeffi- 
cients is similar to boson stars, see for example ||27||. 



in a similarf albeit more involved) form as the metric co- 
efficients (psj) . The same could also be done for the other 
quantities . 
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FIG. 1. Metric coefficients grr — 1 and —gtt — 1 and their 
respective two first terms in expansion, (?,.ro ~ 1, <7ri-2 (top) 
and —gttQ — 1, —gtt2 (bottom), as explained in the text for a 
central value a{0) = 0.2/\/2. The mass of this 0-node oscilla- 
ton is M = 0.5295 mp;/m (see eq. (|27|) and Fig. |), and the 
Schwarzschild metric coefficients are also shown for a solution 
of the same mass. The calculated fundamental frequency is 
Q. = 0.9183 (see eq. (H) and Fig. |). 

The metric functions t'o, i^i, /J.o, Mi are shown in Fig. 
^. It can be verified that the required boundary condi- 
tions are satisfied. As we said before, the central val- 
ues are of the same order than cr(0) while the quantity 
{vi + ^ii) is at least an order of magnitude smaller. These 
are typical behaviors. 

The energy density for the oscillaton (H) is written 



1 



p$(t, a;) — — rn^pim?' 
in 



("^^''V^ sin^ (cut) 
-f e-^'+'^V cos^ (tot) + cos^ {ujt)] (26) 



with as in (Ea). Its values for times ujt — 0,7r/2 
are shown in FigT^. It is easily noticed that /5$(0,a;) — 
p<i>{TT, x), and then the energy density may be expanded 




FIG. 2. Metric functions i/o, vi, fJ-o and /ii for the solution 
given in Fig. nl. 
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FIG. 3. The energy density function p,s>{ijjt,x) (see eq. M) 

for the oscillaton in Figs. Mti (top), and the mass observed at 

infinity (see eq. for different central values (t(0) (bottom). 



Since the metric is asymptotically flat and static, in 
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the sense that the oscillations of the metric are spatially 
confined (see Fig. |l|), the mass seen by an observer at 
infinity may be calculated by the formula p2[ 
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(27) 



The values found using this formula coincides with the 
ADM mass of oscillatons, but the former converges more 
rapidly and is more convenient from the numerical point 
of view. 

The masses of different oscillatons are shown in Fig. ||. 
We can see that there is a maximum mass M^ax — 
0.5522m|,;/m with (Tc(0) ~ 0.235. This critical value 
is important, since it was reported that oscillatons with 
lower central values than (Tc(0) are stable [|. The mass 
calculated in (|2^) is a constant and this implies that the 
masses observed at infinity are the same for all times. 
In fact, we see in Fig. ^ that the oscillaton matches the 
Schwarzschild solution with the same mass|^. 

As we mentioned before, the fundamental frequency 
is given by the asymptotic value exp {v — /i)(oo) = f2~^. 
Recalling that /io(oo) = —vq{oo) ^ and taking into ac- 
count the rapid convergence of function vq (see Fig. |^), 
we obtain the useful formula 



(28) 



The fundamental frequencies of oscillatons are given 
in Fig. ^. Observe that more massive oscillatons oscillate 
with smaller frequencies. 
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FIG. 4. The fundamental frequencies for the same cases 
shown in Fig. |^. 

To finish this section, some lines about the approxi- 
mation (p^. The neglected terms in eqs. (|l|]22|) are 
proportional to I^{xaa') and /3(2/ii), then the approxi- 
mation is good enough if (|a;crcr'|, 2|/ii|) < 1, or equiva- 
lently, cr(0) < 0.5. In this way, the differential equations 



contain only the leading terms of each expansion. But 
this is not exactly true for the KG eq. (p3|), because 
a term proportional to /o(2^i) cos(3wt) ~ C(l) was ne- 
glected. This term is related to the damping term in eq. 
( |rT| ) mentioned above. Thus, some physical information 
of the system is lost when taking just eq. (|l3|), since we 
should have included a time-dependent damping term. 

An easy manner to address this problem is to take more 
terms in the expansion of the scalar field ([l^), and then 
our ignorance on the damping term is transferred to the 
3wt-term. In doing so, the same situation will appear 
with the latter and we will need a 5wi-term. In prin- 
ciple, the procedure have to be continued ad infinitum. 
However, the Fourier series converges rapidly Pj2g|] since 
the damping information transferred to the higher order 
terms becomes smaller in each step. 

The effect of this damping term on the stability of os- 
cillatons can not be treated with the methods presented 
in this paper. But, we should mention that, using a full 
relativistic dynamical evolution of Einstein equations, it 
has been proved that oscillatons are stable configurations 



III. THE STATIONARY LIMIT PROCEDURE 



From the last paragraph, we can say that eqs. (gO[|2J) 
are a first and good approximation for solving oscilla- 
tons. Furthermore, the equations are simple enough that 
we can take a closer look at them and find in which form 
oscillatons and boson stars may be related. 

If we take the weak field limit cr(G) <C 1 (in consequence 
also ^ 1), eqs. (|20|-p3|) may be written as 



Mo = - [1 + - 1)1 , 
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where we have considered that Io{z) ~ 0{1) and Ii{z) ~ 
0{z/2) for z <C 1. The time-independent parts of the os- 
cillaton, eqs. (^jH) separates from the time-dependent 
ones, eqs. ( p4p2 ). 

Eqs. (|29|-|32D can be read directly from eqs. (P-pl]). 
Suppose that the scalar potential and its first derivative 
can be expanded in a Fourier series by taking eq. ( p^ 



y($) = Vb(cr) + V2{(t) cos {2ujt) ■ 
y'($) = Vi[a) COS {ujt) + ... 



(33) 
(34) 



^As it has to be according to Birkhoff 's theorem 
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This is possible if the scalar potential is even. Then, us- 
ing again the dimensionless variables (^8|), eqs. (^9[]3^ ) 
are just 



[x^KoVo{a)/m^ - l] } , 



(35) ^ 
(36) 
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while eq. (^) remains the same. 

The lowest order of approximation to the final result 
of the oscillaton, that in which the time-dependent and 
time-independent parts are separated, can be written us- 
ing a harmonically time-dependent real scalar field and 
expanding the scalar potential in a Fourier series. I shall 
call this as the stationary limit procedure, in order to 
distinguish it from the treatment of section II. However, 
it must be stressed that this procedure will only work for 
potentials with a minimum. Otherwise, the assumption 
of coherent scalar oscillation does not make any sense. 

In Fig. we show the numerical results for the met- 
ric coefficient grr — 1 for an oscillaton using eqs. (Eftl24) 
and for the boson star-like solutions to eqs. ( p^ , |29|-plD, 
with the same central values cr(0) = 0.2/72,0.02/^2. 
The adjusted values of /xo(0), /ii(0) for both cases nearly 
coincide, being almost indistinguishable for the weaker 
field. 

A useful point about the stationary limit procedure is 
that, for the case of a quadratic scalar potential, eqs. ( p9[ - 
|3l| ) are of the same form than those worked for complex 
scalar field boson stars 22,2^]. Indeed, the central 
value (t(0) is also the central value for the equivalent bo- 
son star. 




FIG. 5. Comparison of the numerical results for the met- 
ric coefficient grr ~ 1 for a quadratic oscillaton obtained with 
the method of sec. II (solid lines) and with the stationary 
limit (boson star-like) procedure (dashed lines), for the cen- 
tral values a(0) = 0.2/^ (top), 0.02/^2 (bottom). The 
curves represent g^ro ~ l,<?rr2 as in Fig. 0. 

A complete solution for oscillatons should have more 
Fourier terms in (n3[), that is, we should add more field 
contributions ||5|,p9| . In general, the curves for an oscil- 
laton are (a little bit) smaller than the those of a boson 
star represented in this case by the stationary limit pro- 
cedure. This is a general fact since, unlike boson stars, 
oscillatons have an extra damping term in the KG equa- 
tion ([III). Then, it would be found that oscillatons are 
in between boson stars (stationary limit procedure) and 
the results found in the previous section. 

However, the boson star-like equations are easier than 
a full treatment in Fourier series. Thus, the stationary 
limit procedure can be used to determine an upper limit 
for oscillatons and to provide, at least, order of magni- 
tude results of interesting quantities. The results become 
exact when dealing with weaker fields. From now on, 
these are the general criterions we will take when working 
with the stationary limit procedure (see next section). 

For instance, the maximum mass for a (O-node) boson 
star Mbs c = 0.633mp;/m is reached when a{0) = 0.268 
p2| , while the maximum mass for a (O-node) oscillaton 
is A/$,c = 0.5522m,%i/m with cr(0) ~ 0.235, see Fig. |. 
As we said before, the value of the maximum mass given 
by the boson star-like equations is an upper limit and of 
the same order of magnitude than the exact value . 

An interesting point now arises. From Fig. |l| we 
observe that the metric for an oscillaton is static from 
a; = 10 up. From Fig. ^, we also observe that it is dif- 
ficult to distinguish between an oscillaton and a boson 
star far from the center, and in the weak field limit, the 
results are almost the same. This result suggests that 
the oscillaton is some kind of boson star plus an intrinsic 
oscillation (albeit this is not exact for non-polynomial po- 
tentials, as we will see in section V). The two objects are 
different only for an observer near the center. This last 
fact would imply that oscillatons could compete with bo- 
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son stars as astrophysical objects of interest, since some 
useful properties of boson stars arises in situations far 
away from the boson star |^^, but the picture becomes 
more interesting if we move closer to the origin. 

Concerning the possible Newtonian limit for oscilla- 
tons, we can see that such limit can not be obtained as 
a post-Newtonian expansion, as it was first pointed out 
in 1^. The stationary limit procedure seems to suggest 
that the static parts of the metric could have a kind of 
Newtonian limit as boson stars ||2l|,|2^,^ , but the intrin- 
sic oscillation, the very imprint of oscillatons, will not 
disappear |^. 



IV. SELF-INTERACTING OSCILLATONS 

Another class of interesting oscillatons are those with 
self- interacting scalar fields. Since these cases are more 
involved than the quadratic one, we will consider them 
within the stationary limit procedure only. 

Thus, the aim of this section is to obtain the qualita- 
tive behavior and approximate (at least in the order of 
magnitude) solutions to the leading terms in the Fourier 
expansions for the oscillatons. As it was said before, the 
true solutions must be smaller^ but of the same order of 
magnitude. 

As we will see below, self-interacting scalar fields pro- 
vides with a wider range of properties that could be also 
of astrophysical interest. 



A. The quartic case: V = (m^/2)$^ -|- 
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(41) 



Mo 



X 



{ 1 + e^^+i"" [x^ (0-2 + (A/2)a4) - l] } , (42) 
a' - ^i'^ + (cr + Acr^) - e^'^V (43) 



/I'l = e'^o+^V [xa + (2/3)Aa;(T3 - cr'l 



(44) 



The parameter ( |39D makes the last equations corre- 
spond to the case of boson stars [pOl 



25 1 , and then the sim- 



ilarity between oscillatons and boson stars is preserved at 
this coupling. Following a similar argument as in ]20[ | , the 
self- interaction may be ignored if A <C 1. But if A ^ 1, 
the oscillatons obtained may differ significantly from the 
non-interacting ones. 



0.045 

0.04 
0.035 

0.03 
0.025 

0.02 
0.015 

0.01 
0.005 






A=0 




A=100 




A=204.6 







10 



20 



30 



60 



70 



This kind of potential has been considered as a model 
of bosonic dark matter [|[0J§1 . In this case, we will also 
set the dimensionless parameter 



A: 



3g 



(39) 



with g/A the coefficient of the quartic term. Note that 
this parameter is different than that used for boson stars 
pO| , but it is still a rough estimation of the ratio between 
the quartic and quadratic interactions in the potential. 

A complete treatment for this case is beyond the pur- 
pose of this paper. However, the stationary limit pro- 
cedure can still give us enough information about the 
leadingterms of the oscillaton functions. Taking this 
limit (pq-pq), the equations to be solved are 



vi = xaa — /ii. 



(40) 




FIG. 



10 20 30 40 50 60 70 80 90 
X 

Comparison of the results for the metric coefficients 



OrrQ — i (top), grr2 (bottom), for the case of an oscillaton with 
different quartic self-interactions within the stationary limit 
procedure. The central value is ct(0) = 0.02/i/2. 



■^The extra damping term in the KG equation ( p4| ) does not 
depend on the strength of the scalar potential, then the results 
must not be much smaller. 
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Typical numerical results of an oscillaton with a quar- 
tic self-interaction within the stationary limit procedure 
are shown in Fig. As in the quadratic case, we 

can calculate the mass and frequency of those solutions: 
M = {0.736, 0.4225} TO|,;/m and Q = {0.9787, 0.9857} 
for A = {204.6, 100}, respectively. These numbers can 
be compared with the <i>^ case, M = 0.235 mp;/m, fl = 
0.9904. The self-interacting oscillatons can be more mas- 
sive and oscillate stronger but with a smaller frequency. 

We can take advantage of the similarity with boson 
stars to approximately determine an upper limit for the 
maximum mass of an oscillaton with large-A. The result 
would be given by pOt] 



Mrr^ax ^ 0.2\/Am^i/ m 



(45) 



where A is given by (p9|). 

Another interesting point has to do with the critical 
value (Tc(0) at which an oscillaton reaches the maximum 
mass. It is also important because oscillatons would be 
stable only for values cr(0) < <Tc{0) as it is the case for 
boson stars | |25[ |. 

Recalling the case of boson stars, the critical value 
becomes smaller for larger A. For instance, (Tc(0) — 
{0.1,0.07} for A = {100,204.6} respectively, see Fig. 1 
in |2^]. Moreover, for large-A, the critical value is ap- 
proximately given by 



(46) 



as suggested by Fig. 4 in pO| |. 

This result can be taken as an upper limit of adO) for 
oscillatons. That is, we can say that there must be sta- 
ble oscillatons only for central values smaller then (46). 
In both cases, boson stars and oscillatons, we need 
weaker fields for supporting stable self-interacting con- 
figurations. 

Other useful information may be obtained from boson 
stars |2^,^ and then approximately applied to oscilla- 
tons through the stationary limit procedure. In a similar 
way as in boson stars, quartic-oscillatons would provide a 
new range of configurations with more interesting prop- 
erties than the quadratic ones. 



B. Non-polynomial potentials: 

!/($) = Vo [cosh (Av^$) - 1] 

In a series of recent papers, our group have been inves- 
tigated the hypothesis of scalar field dark matter in the 
Universe [pl]-p9[. It has been shown that a minimally 
coupled scalar field <I> endowed with a scalar potential 



y($) = Vo [cosh {\y^<P) - 1] (47) 

mimics quite well all known and desirable properties of 
cold dark matter at cosmological scales 
over, the free parameters of the potentia' 



TlJTI]. More- 
Vb, A can be 



fixed from cosmological observations jlj] . The key point 
for this model is that the scalar potential possesses a min- 
imum and that the scalar potential behaves as ~ <I>^ at 
late times in the evolution of the Universe. This scalar 
potential belongs to the family of scalar potentials of the 
form 



U($) = VbMnh(aV^$)] 







which has been studied in quintessence ||lC 
[lO[|l|-|l|,0,@. 



(48) 

511 and dark 



matter 

What kind of oscillatons would appear with ([iTf)?. 
First of all, the self-interaction is a very important issue 
JlSf . At first sight, we may take (47) with just a quartic 
coupling as in the previous section, with g — kqui^X^ /3l 
and the scalar field mass defined by — kqVoX^- From 
eq. (^9|), we have then that 



2 



(49) 



For instance, taking a scalar dark matter model with a 
cosh potential (l^Jl^ cosmological observations suggest 
A — 20.28 (A ~ 205.6). In this particular case, we would 
have strong self-interacting oscillatons like those shown 
in Fig. |. 

The full treatment of potential ( p7| ) is even more com- 
plicated than the quartic one. In order to have approxi- 
mate information about the leading terms of the solution, 
let us take the stationary limit procedure ( |35[]3^ ). Con- 
sidering that 

[cosh (2Act cos (cjt)) - 1] = (/o(2A(t) - 1) 

+2/2(2Act) cos (2ujt) + 
sinh (2Acr cos (ujt)) = 2/i(2Acr) cos (ujt) + 



the Einstein-Klein-Gordon equations read 



— xaa — /ii. 



in 

X 



-(/o(2Aa)-l)-l 



(50) 
(51) 

(52) 



1 - - Mo 



-e''»+^»/i(2ACT) -e^'^'cr (53) 
A 



Ml 



(54) 



It is worth to notice that the mass of the boson also 
sets the scale of the solution as in the cases before, albeit 
a more natural choice for this case could be the quantity 
VkqVo ||1^. The appearance of Bessel functions when 
expanding the scalar potential also shows that all cou- 
plings continue participating. Unfortunately, the case 
of a cosh-boson star has been treated before |2^ and 
then we see that the similarity among cosh-oscillatons 
and cosh-boson stars does not exist anymore. 

Since we are mostly interested in stable cosh- 
oscillatons, we should first estimate CTc(O). From the 
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<I>^-oscillatons we know that the stronger the interaction, 
the smaher the critical value adO). Then, the value (^ ) 
must be an upper limit for the cosh case. Therefore, for 
cosh-oscillatons we obtain 



crc(O) ~ V2X' 



(55) 



But, observe that for these values (and smaller ones) 
the cosh potential (j^ can be treated perturbatively. 
That is, eqs. ( ^o[]5^ ) can be worked in the Hmit of eqs. 
( ^Ofji^ ). Therefore, as long as we consider only stable 
oscillatons, the cosh potential gives (almost exactly) the 
same results as the quartic potential does. 

Going further, the maximum mass for the cosh oscilla- 
ton must be of order M^ax — 0.1Amp;/m (see eq. (^5|)). 
This formula can be written in a more convenient form 
by taking the definition of the scalar field mass. Then, 
the maximum mass for a cosh oscillaton is approximately 
given by 



0.1 



VkqVo ' 



(56) 



which in turn means a A-independent value. 

The validity of the last formula and the existence 
of stable cosh-oscillatons have been recently confirmed 
using a full numerical evolution of the Einstein-Klein- 
Gordon equations p^. 



V. CONCLUSIONS 

The simplest approximation for solving the coupled 
Einstein-Klein-Gordon equations for a spherically sym- 
metric -oscillaton were presented, by taking an har- 
monic time-dependent scalar field. Non singular and 
asymptotically flat solutions are obtained taking a 
Fourier expansions of the differential equations. These 
solutions are also asymptotically static, so the metric os- 
cillations are confined to the region of the oscillaton. Us- 
ing this last fact, we could determine the mass and funda- 
mental frequency of some configurations. The questions 
concerning the stability and formation of these oscilla- 
tons were answered based in previous publications. 

We presented a method called the stationary limit pro- 
cedure, which gives the lowest order of approximation 
to the leading terms in the Fourier expansions. In this 
limit, the differential equations for the oscillaton are bo- 
son star-like and the similarity between these objects be- 
comes manifest. Taking this similarity, we could deter- 
mine an upper value and an order of magnitude estima- 
tion for some critical parameters of oscillatons by com- 
paring them with boson stars. 

The cases of quartic and cosh self interactions need a 
more complicated numerical treatment. However, they 
were studied in the stationary limit procedure. For the 
quartic case, boson star-like differential equations arise 



again. Thus, the solutions given in the literature for bo- 
son stars apply for oscillatons, at least at the level of dif- 
ferential equations. In the non-polynomial case of a cosh 
scalar potential, the equations are not boson star-like, 
but the results are quite similar to the case of quartic 
self-interaction, as long as we refer only to stable con- 
figurations. About the stability of these oscillatons, we 
determined the central value adO) at which the maxi- 
mum mass is reached. In the quadratic case, there is 
evidence that the oscillaton is unstable for central val- 
ues larger than (Tc(0), and we assume the same for the 
self-interacting ones. Whether this really happens is still 
matter of current research and the results will be pub- 
lished elsewhere lT^p9[ |. 

There are some issues to be investigated. First of all, 
it is desirable to have a more useful and complete nu- 
merical treatment to find the oscillaton configurations, 
one in which the non-linearities of the EKG equations 
could be minimized as much as possible j2^. A simi- 
lar step should be made to include self-interacting os- 
cillatons. This would allow us to exactly determine the 
reliability of the stationary limit procedure. 

As it was shown through the stationary limit proce- 
dure, it must be possible to find a kind of Newtonian-like 
limit of oscillatons. Apart from the simplicity and the 
appealing scaling properties that may appear (for boson 
stars and their Newtonian limit see |2l| , |23| , ^ , p2| ), such 
limit would be also of cosmological and galactic inter- 
est, like in the case of Newtonian boson stars ||2^ , ^ , ^ . 
Work in that direction is in progress |^^. 

Finally, scalar dark matter would provide of very in- 
teresting gravitationally bound objects like oscillatons. 
These relativistic systems, with no Newtonian counter- 
part, would appear from the gravitational collapse of a 
real scalar field, and their properties may be directly con- 
nected to cosmological parameters 13 -19 . This addi- 



tional reason makes oscillatons worth investigating. 
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